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The dynamics and internal structure of doping fronts in organic semiconductors are investigated 
theoretically using an extended drift-diffusion model for ions, electrons and holes. The model also 
involves the injection barriers for electrons and holes in the partially doped regions in the form 
of the Nernst equation, together with a strong dependence of the electron and hole mobility on 
concentrations. Closed expressions for the front velocities and the ion concentrations in the doped 
regions are obtained. The analytical theory is employed to describe the acceleration of the p- and 
n-fronts towards each other. The analytical results show very good agreement with the experimental 
data. Furthermore, it is shown that the internal structure of the doping fronts is determined by the 
diffusion and mobility processes. The asymptotic behavior of the concentrations and the electric 
field is studied analytically inside the doping fronts. The numerical solution for the front structure 
confirms the most important predictions of the analytical theory: a sharp head of the front in the 
undoped region, a smooth relaxation tail in the doped region, and a plateau at the critical point of 
transition from doped to undoped regions. 



I. INTRODUCTION 

Organic semiconductors (OSCs) demonstrate a num- 
ber of interesting properties, which distinguish them 
from crystalline inorganic semiconductor*^. Among the 
most intriguing features is the possibility of electrochemi- 
cal doping by means of reversibl e red ox reactions, demon- 
strated in a number of materials^. The electrochemical 
doping transforms the OSCs from essentially insulating 
state to a metallic-like state with the conductivity in- 
creased by many orders of magnitude^. This transfor- 
mation is accompanied by a considerable change of the 
material properties including color, p hotolu minescence 
capability, volume, and surface energy^Hl. The elec- 
trochemical doping can be performed straightforward in 
situ, when an OSC coated on a metal electrode is in con- 
tact with an electrolyte and is subjected to an appro- 
priate electric potential. In the transformation process, 
electronic charges (electrons and holes) from the elec- 
trode are injected into the OSC and subsequently elec- 
trostatically compensated by an influx of respective ions 
from the electrolyte (cations and anions 

The opportunity for tuning the electronic and opti- 
cal properties of OSCs has triggered a number of stud- 
ies at the fundamental as well as practical aspects of 
doping. These works in turn paved the way to nu- 
merous applications of this p rocess in novel electronic 
and photonic devices^^^HXSI a classical example of 
the electronic devices utilizing the electrochemical dop- 
ing is polyme r based light-emitting electrochemical cells 
(LECsJPESEElHIID A LEC comprises an organic semicon- 
ductor in a form of conjugated fluorescent polymer, which 
is blended with a solid state electrolyte. The blend of 
the conjugated polymer and the electrolyte forms the 
active material of the LEC and is sandwiched between 
two electrodes. When an electric potential applied be- 
tween the electrodes exceeds the band gap potential of 
the conjugated polymer, then the doping transforma- 



tion starts with injected holes and electrons forming 
the p- and n-type doped regions close to the respec- 
tive electrodes^. The doping process in planar LECs 
can be visualized under ultra-violet illumination since 
doped OSCs display very high rate of photoluminescence 
quenching. The doped material in LEC is observed as 
dark regions quite distinct from the original undoped 
substance. By employing the ultra-violet visualization, 
it was demonstrated experimentally that doping trans- 
formation in OSC develops in the form of two localized 
fronts of p- and n-type doping, which emerge at the elec- 
trodes and propagate towards each othef^Hli] When the 
fronts meet, the two doped regions form a p-n junction, 
which emits a visible light. The purpose of the present 
work is to provide a theoretical model for the front prop- 
agation in LECs prior to the development of the p-n junc- 
tion. 

A number of interesting theoretical works has been 
devote d to charge dynamics in LECs and other OSC 
devices^EMUl i n particular, the authors of Refs. 1221 
and 1251 investigated a stationary light-emitting p-n junc- 
tion as a final state of the doping process in LECs. At 
this stage of the process, the whole OSC is already con- 
verted to the state with high conductivity. However, the 
non-steady problem of doping front dynamics and struc- 
ture prior to the formation of the p-n junction is much 
more difficult to analyze, since it involves transition of 
the OCS from the undoped weakly conducting state to 
the metallic- like doped state within the front. As we 
show in the present paper, the problem includes not only 
the electrodynamic issues of the OSC plasma motion, 
but also some questions of thermodynamics and quan- 
tum mechanics, which are still waiting for an answer. 
The problem of front dynamics was addressed in [20] and 
[211 from the empirical point of view of total OSC con- 
duction: the purpose was to analyze the experimental 
data without investigating the complicated internal front 
properties. At the same time, considerable progress has 
been achieved in the adjacent problem of front dynamics 
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in one-electrode OSC de vices, su ch as e.g. electrochemi- 
cal sensors and actuator d 24 * 26 * 27 ! For example, the recent 
work by Wang et alp3 provided discussion of all elements 
required for front formation in the one-electrode devices. 
The model proposed in Ref. [37J included the Poisson 
equation, diffusion and mobility (migration) of holes and 
ions (cathions), taking into account possible nonlinear 
(non-Fickian) character of the transport coefficients. The 
Nernst equation for light particles, holes, in the region of 
high conductivity in that case follows from the diffusion- 
mobility model under the condition of quasi-equilibrium. 
In a sense, the process studied in Ref. [57] is just the op- 
posite to the doping fronts propagating in LECs. In the 
doping process in LECs, holes and electrons are injected 
into the active material by an externally applied electrical 
field, which helps them overcoming a certain thermody- 
namic barrier, and dope the active material. The ions, on 
the other hand, give way to the light charges, and com- 
pensate the excessive charge thus avoiding generation of 
strong internal electric fields. In the devices presented by 
Wang at alp2 charge motion goes in the opposite direc- 
tion: the cathions advance together with the front, while 
holes retreat leaving the material. With some caution, 
the fronts of electrochemical transformation in LECs and 
in the one-electrode devices^ may be compared to uphill 
climb and downhill glide, respectively. Still, the physi- 
cal understanding of front dynamics in the one-electrode 
devices as presented irP may be helpful in constructing 
models of the doping front propagation in LECs. In par- 
ticular, recent paper 2 ^ demonstrated that the theoretical 
description of doping fronts in LECs requires not only the 
common diffusion-mobility set of equations, but also non- 
linear concentration-dependent transport coefficients for 
holes/electrons and the thermodynamic injection barrier 
for the light charges in the form of the Nernst potential. 

The present paper develops the ideas of Refs! 27 - 12 ^. Us- 
ing the theoretical model proposed irPll, we study dynam- 
ics and internal structure of the doping fronts in LECs. 
We derive compact analytical formulas for the front ve- 
locities depending on the electric field together with the 
ion concentrations in the doped regions. On the basis 
of the analytical theory we describe acceleration of the 
p- and n-fronts approaching each other. The analytical 
results show very good agreement with the experimental 
data. We show that the internal structure of the dop- 
ing fronts is determined by diffusion and mobility pro- 
cesses. We study analytically the asymptotic behavior 
of the concentrations and the electrical field inside the 
doping fronts. We also solve numerically for the front 
structure, and confirm the most important predictions of 
the analytical theory: a sharp head of the front in the 
undoped region, a smooth relaxation tail in the doped 
region, and a plateau at the critical point of transition 
from doped to undoped regions. 

The paper is organized as follows: In Sec. II we in- 
troduce the basic equations of the model. In Sec. Ill we 
study properties of the doping fronts considered as sur- 
faces of discontinuity; we find the front velocity and ion 



concentrations in the doped regions behind the fronts and 
then describe front acceleration in LECs. Section IV is 
devoted to the internal structure of the fronts. We discuss 
characteristic length scales of the process, the condition 
of quasi-neutrality, and the necessity of an injection bar- 
rier for the system of equations. We investigate asymp- 
totic analytical solutions to the equations in the specific 
zones of the front and, finally, we numerically solve the 
whole set of equations. The paper is concluded by a brief 
summary in Sec. V. 



II. THE MOBILITY-DIFFUSION APPROACH 

In general, the dynamics of the ions is determined by 
equations of force balance 

nm — = —qnSI(j) — ksTVn nrav, (1) 

dt t 

where q = ±e is charge of positive/negative ions, m is ion 
mass, n is concentration, <f> is potential of an electric field 
with E = — V0, kg is the Boltzmann constant and T is 
the polymer temperature, which may be taken constant 
and uniform. The last term in Eq. ([I]) takes the aver- 
age contribution of collisions into account. In the case of 
organic polymers, the collisions dominate over the iner- 
tia terms, which allows the mobility-diffusion approach, 
giving the velocity according to 



n± 



(2) 



where the mobility is given by fi — er/m, labels ± cor- 
respond to positive and negative ions, respectively, and 
the diffusion coefficient D is related to mobility using 
the Einstein relation fi = eD/ksT. When substituting 
velocity from Eq. Q to the continuity equations 



dn± 
~dt~ 



•V-(n±v±) = 0, 



(3) 



we arrive to the mobility-diffusion model for ion motion 
according to 



dn± 

Ik 



V • [±fi ± n±V</> + D±Vn±] = 0. 



(4) 



The equations for electrons and holes may be presented 
in a similar way, though in the case of light charges one 
has to take into account an injection barrier </>jv for the 
electrons and holes in the transition from the doped re- 
gions to the undoped ones 



dn h 
dt 



dn e 
~dt 



V • [nh.n tl V( 



V ■ [-/x e n e V( 



) + D h Vn h ] = 0, (5) 



D e Vn e ] = 0, (6) 
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p-doping front propagating with velocity U p along the x- 
axis in a static and uniform electric Eq = const, created 
in the undoped region by external sources. Looking for 
the p- front solution in the form = ^(x — Ut) we reduce 
Eqs. Q, @, Q to 



_d_ 

dx 



-n h U p + n h n h (E - E N ) - D, 



d_ 

dx 



-ri—Up — ri-fi-E — D. 



dn h 
dx 



= 0, (9) 



dri- 
dx 



= 0, 



(10) 



FIG. 1. Schematic of the doping process. 



d_ 

dx 



-n + U p + n+fi+E — D 



dn + 
dx 



= 0, 



(11) 



where the labels "h" and "e" stand for holes and elec- 
trons, respectively. The barrier originates from quan- 
tum and thermodynamic effects, and the interpretation 
of (ft 7v in terms of an electrostatic force should therefore 
be done with caution. We discuss the necessity of <Pn for 
self-consistent description of the doping front structure 
and the particular form of this term in the Section IV. 
Meanwhile we point out that this term is non-zero only 
inside the doping fronts, and it should turn to zero in 
the uniform undoped and doped material ahead of the 
fronts and behind the fronts, respectively. The electric 
field obeys the Poisson equation 



V 2 </> = (n_ + n e - n + - n h )e/e Q . 



(7) 



Still, in Section IV we demonstrate that the condition of 
quasi-neutrality 



n + — rih = 



(8) 



holds with a very good accuracy in LECs. We stress 
that condition ^ does not mean constant electric field 
everywhere. Instead, it implies that even tiny local devia- 
tions from zero net charge lead to extremely large electric 
fields. 



III. 



PROPERTIES OF DISCONTINUOUS 
DOPING FRONTS 



It was demonstrated experimentally^ 20 * 21 ! that the 
doping process in OSCs happens in a form of two fronts 
propagating towards each other as shown schematically 
in Fig. 1. A p-doping front populates the semiconductor 
with holes, while the n-front makes the semiconductor 
rich with electrons. Both fronts propagate due to exter- 
nal electric field created by the potential difference be- 
tween the ends of the semiconductor film. In the present 
section we consider planar p- and n-fronts as propagat- 
ing discontinuity surfaces, which transform the original 
undoped semiconductor to a doped conducting material. 
To be particular, we start our analysis with a stationary 



dE 
dx 



= e(n + - n_ + n h )/e 0l 



(12) 



where Ejv = — V0at is the effective electric field related 
to the quantum-thermodynamic injection barrier. Still, 
in the limit of discontinuous doping fronts considered in 
the present section, the barrier term does not play any 
role, since it is zero in the uniform regions both ahead 
of the front and behind it. Equations ^ - ( fl2] ) may 
be integrated analytically across the front. We designate 
values ahead of the front by the label "0" and values be- 
hind the front by "1". We also take into account that 
the hole concentration is zero ahead of the front and it 
reaches some known finite value n\h behind the front. 
The initial concentrations of positive and negative ions 
are equal and known, no- = "0+ — n a- Then the pa- 
rameters in the doped and undoped regions are related 
by the integrals of Eqs. ([9|-( 11 ) according to 



-U p + fihiEi = 0, 



- n 1+ U p - n 1+ [i + Ei = -n Q U p - n Q fi + E Q , 



nx-Up + nx-[i-E\ = —n U p + n ^i^E , 



(13) 



(14) 



(15) 



Equation ([13j) specifies the electric field in the doped 
region as E\ = U p /^hi, which may be neglected tak- 
ing into account high mobility of the holes in compari- 
son to the ions in the doped material: E\ = U p /fihi 
{n±/ Hhi)Eo -C Eq. Equation (12 1 in the uniform doped 
region determines zero net charge as 7ii_ — n 1+ = n\h- 
Then equations (14), (15) specify the front velocity and 



ion concentration in the doped region as 



Up = — (fi+ + fi-)E Q , 
nih 



ni- = n +ni h 



n\+ = n - n lh 



/'- 



fl + + [A— 

M+ 
/i+ + H- 



(16) 
(17) 
(18) 
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Similar formulas may be obtained for the n-doping front, 
i.e. 



U n = — + (i-)E , 
n 2e 



n 2+ =n Q + n 2e - 



n 2 _ = n - n 2c 



(19) 



(20) 



(21) 



where the label "2" designates the conducting substance 
behind the n-doping front. We note that the p- and n- 
type doping fronts propagate with different velocities re- 
lated by 



U n = —U p = (3U P . 



(22) 



Taking data for electron and hole concentration n\h = 
8.6 • 10 25 m ~ 3 , nie = 1.3 • 10 26 m- 3 obtained in the 
experiment j 16 l 21 l2£l we find the velocity of the n-front to 
be less than that of the p-front with (3 — 0.661. 

Next we consider a semiconductor of finite size with 
a voltage <fio applied to the end electrodes and two pla- 
nar doping fronts moving towards each other. The initial 
distance between the fronts is Lq. Since the distance be- 
tween the fronts L(t) decreases in time and the potential 
difference is fixed, then electric field Eo(t) in the gap be- 
tween the fronts grows in time and the fronts accelerate 
according to Eqs. (16), (19 1. We designate position s of 

In experiment d 16 | 21 | 2 SI, the 



t) 



the fronts by X p (t), X„ 
n-front (electron doping) starts later, after a time delay 
t e ; we thus have two time intervals in the solution: t < t e 
and t > t e . We start with the first interval, t < t R , when 
only the p-front propagates. In that case Eq. 
the p-front velocity is reduced to 



(16) for 



dX p 
~dt~ 



n 
nih 



,9o _ 
L 



U, 



pO 



1 — X p / Lq 



(23) 



where U p q is the initial velocity of the doping front, 
tegrating ( 23 1 we find 



In- 



X,, 



E = 1- Jl-2U p0 t/L . 



(24) 



At the moment t e , when the n-front starts, the p-front is 
already at the position X p (t e )/L n = 1 — — 2U p ot e /L , 
and we have the reduced distance between the fronts 
L e = L(t e ) = L Q — X p {t e ). Next, we consider the in- 
terval t > t e , when both fronts propagate. In that case 
the distance between the fronts varies as L = X n — X p 
with 



dX. 



U, 



pO 



dt 



L/L 



dXn 

dt 



po 



L/Lq' 



so that 



dL 
~dt 



-(1 



u 



pO 



L/L - 



(25) 



(26) 
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FIG. 2. Positions of the doping fronts versus time as predicted 
by the analytical theory (solid lines) and obtained experimen- 
tally (markers). The uncertainty bars indicate the difference 
between the fastest and the slowest parts in the multidimen- 
sional front structure. 



Then, for t > t e , we obtain the solution 



L(t) = yj 'l\- 2(1 + /3)Up L (t-t e ), (27) 



Xp Xp (t e 



1 + P 



[Le-L(t)}, 



1 + /3 



[L e -L(t)]. 



(28) 



(29) 



Figure 2 shows a comparison of the positions of the pla- 
nar p- and n-fronts, as found analytically through Eqs. 



(28), (29), to the experimental data. The markers show 



the average position of the fronts; the error bars indicate 
the difference between the fastest and the slowest parts 
of the two-dimensional (2D) front structure (width of the 
2D front brush) . Such 2D structures develop due to front 
instabilities, which are beyond the scope of the present 
paper and which will be presented elsewhere^. Figure 
2 indicates that the analytical ID solution describes the 
dynamics of the backside of the front structure quite well 
(within the 10% accuracy of the experimental data). At 
the same time, the leading parts of the front move no- 
ticeably faster than the ID theory predicts. The last fact 
demonstrates an important role of the 2D instabilities in 
the doping front dynamics, which is the subject for the 
future work. 
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IV. INTERNAL STRUCTURE OF A DOPING 
FRONT 

A. Basic equations 

In the present section we investigate the internal struc- 
ture of a stationary doping front driven by a constant 
uniform external electric field Eq. We start by study- 
ing the p-front, since analysis of the n-front is similar. 
The internal structure of the doping front is controlled 
by diffusion in Eqs. ([9])-( 11 1 integrated within the front 



n h U p + n h n h {E - E N ) = D, 



dx 



(tiq — n-)Up — /j_(?i_i? — uqEq) = D. 



(n - n + )U p + u+(n + E - n E ) = D. 



dri- 
dx 

dn + 
dx 



(30) 



(31) 



(32) 



A characteristic length scale L p involved in the problem 
is determined by the slow diffusion of ions, thus we may 
choose, e.g., D_ to define 



Ln 



D-n lh 



U p (fj,+ + U-)n E ' 



(33) 



m 

,2 A/ j> 



of 

3 



As an exa mple, taking the experimental data 
RefeGSEEHfor Ulh = 8 . 2 . 10 25 m" 3 , n = 3.1 • 10 26 
H + = 1.0 • 10- 10 m 2 /Vs, fj,_ = 2.2 • 10- 10 m 2 /Vs, 
360K, Eq — 3 ■ 10 3 V/m, we can evaluate the characteris- 
tic initial p-front velocity as U p o « 3.6 • 10 _6 m/s and the 
initial length scale as L p q 1.8 • 10 _6 m. Still, one has 
to remember that L p does not portray the full structure 
of the doping front. We show below that there are sev- 
eral different characteristic length scales within the front 
from the undoped to the doped region. 

We may also estimate the typical deviations from 
quasi-neutrality within the front through 



n + — n_ + ri/j e dE £qE q 
uq eno dx eL p ono 



2.9- lO" 10 < 1, 



which allows replacing the Poisson equation Eq. (12 1 by 
the condition of quasi-neutrality 



(34) 



n_ 



(35) 



everywhere within the front with very good accuracy. 
As a comparison, the electrochemical fronts in the am- 
bipolar devices involve noticeable deviations from quasi- 
neutrality as shown in RefP-2 

The next important question concerns the terms speci- 
fying the injection barrier, cf>N, En, in Eqs. ([5| and (30 1. 
Below, we demonstrate that it is indeed impossible to de- 
scribe the doping front structure without such a barrier 
term. Let us consider the doped region at the back side 



of the front with /i^ ^> u±. Without the barrier term, 
one has the force-balance equation 

n h m -ji = -en h \/(f> - k B T\7n h v. (36) 

at t 

for holes similar to Eq. ([IJ. Within the drift-diffusion 
approximation, the inertial term nyjndvjdt in Eq. (36) 
is negligible. The characteristic velocities of the doping 
process are controlled by relatively low ion mobility, Eq. 
Therefore, by order-of-magnitude the collision term 
tjq. ( 36 1 is also negligible in comparison with the 



(16 



electric term in the doped region 
n h m m 

v oc — no(fi+ + /i_JV0 cx 

T T 

^ + — ertfr <C en,hV(j). 



(37) 



Then, in the doped region, Eq. (36) should describe 



hydrostatic equilibrium for the holes, with balance be- 
tween electric and pressure forces proportional to — 
and -Vrifc, respectively. However, in the geometry of 
the doping process, the electric and pressure forces point 
in the same direction (e.g. to the right in Fig. 1) and 
therefore cannot balance each other. Without the bar- 



rier term, Eq. ( 36 1 suggests that holes are freely acceler- 



ated into the semiconductor, which contradicts the very 
essence of the doping process. In order to obtain a doping 
front, one has to balance the electric and pressure forces 
in Eq. (36) by a counter-term, which takes into account 



the thermodynamic barrier of the doping process 



dw 

nhm— = -en,jV( 



) - k B TVn h 



nyjn 



(38) 



dt t 

It has been suggested in RefPSlfo consider </>jv in the form 
of the Nernst potential, given by 



<j) N (n h ,T) 



ln 



nh,oa - n h 



(39) 



where n^.oo is the maximal possible concentration of 
holes. In general, n^oo may be larger than the value 
7i/ji of the experimentally observed hole concentration 
behind the p-doping front. For this reason, in the present 
paper we introduce also a numerical parameter / as 
nhoo = ffihi with / > 1. We will show that qualita- 
tive properties of the doping front do not depend on the 
parameter /. As mentioned above, the Nernst potential 
should be viewed as an injection barrier, and is only de- 
scribed as an electrostatic contribution out of convention. 
Moreover, the important dynamical aspects of the Nernst 
potential enter Eq. (36) in the form of a gradient. Since 



the inclusion of the Nernst potential is formally valid only 
in the highly doped region, it can therefore be set to an 
arbitrary constant value in the undoped region, i.e. 



<j>N(n h ,T) = 0. 



(40) 



A continuous description of the doping front structure 
also requires a continuous change of the Nernst poten- 
tial from Eq. (39) to Eq. (40). Since, unfortunately, 



(> 



there is no theoretical thermodynamic model describing 
such a transition at present, we introduce a dimensionless 
phenomenological function ij) with the injection barrier 
potential given by 



0N(n h ,T) = ip 



In 



nh 



(41) 



where ip = 1 in the doped region and ip — in the un- 
doped region. A more detailed form of the function ip 
will be discussed below. Taking into account the Nernst 
potential, Eq. ( 30 ) is modified in the doped region ac- 
cording to 



n h U p + n h [i h E = D h I 1 - i> 



fn h \ 



fnhi 



dn h 
dx 



(42) 



The set of Eqs. (31), (32 



(35), (42) determine the in- 



ternal structure of the p-doping front. 

Another important feature of the system is a strongly 
nonlinear dependence of the hole mobility on concentra- 
tion Hh{nh) found experimentally^. Because of this de- 
pendence, hole mobility is much larger than the ion mo- 
bility in the doped region [ih 3> fJ>± , but it becomes much 
lower (x h <C fx± in the undoped region where n^/n^i <C 1. 
The nonlinear dependence is consistent with basic under- 
standing of doping as the process of increasing mobility 
of charge carriers. Reference^ makes use of the following 
empirical fit for the hole mobility 



Hh = 3.85 x 10" 



1 + tanh ( 26.6— - 4.3 
n 



A similar property holds also for electrons. 



m 2 /Vs. 
(43) 



B. Dimensionless equations 

In order to simplify the analysis, we introduce dimen- 
sionless variables for the coordinate, concentrations and 
electric field: £ = x/L p , a± = n±/no, ah — Tih/no, 
s = E/Eq. Furthermore, we introduce two parameters 
specifying ratio of ion mobilities and the front velocity 
according to 



5 = fi-/fx + , C = fi-E /Up . 



(44) 



In particular, we have S = 2.2 for the active material 
used in the experiments of RefsP^2U2SI l n the dimen- 



sionless form, the governing equations (42), (31) 
( 35 ) become 



(32), 



da- 
da + 



= 1 + C - a- (1 + Ce) , 



6-C-a+(6~Ce) 



(45) 



(46) 



(47) 



a + — a- 



ah 



0. 



(48) 



where the dimensionless hole mobility and the effective 
dimensionless hole diffusion are 



M h = 175[1 + tanh (26.6^ -4.3)], 



Te// 



= M h 1 



Ipahoo 

ahoo - a h 



(49) 



(50) 



The numerical solution also requires boundary condi- 
tions in the doped region, as specified in Sec. Ill under 
a simplifying (though realistic) assumption of infinitely 
large hole mobility in the doped region, /i/ji 3> fi±. Still, 
it is also useful to determine the boundary conditions 
without employing such an assumption. The hole concen- 
tration in the doped region is known from experiments. 
Thus, we have to find the ion concentrations in the doped 
region and the exact value of the C-parameter in Eq. 
(44), taking into account a finite (though small) electric 
field E\. Since all derivatives in Eqs. |45|)-(47) are zero in 



the doped region, then we obtain the following boundary 
conditions 



£i 



1 

CM lh ' 



0-1- 



l + C 
l + Cei 



an 



6-C 
6-Cei 



(51) 



Substituting (51) into (48) we find the parameter C 



C 



aih 



1 



M u 



1 



M 



a,. 



1 



M- 



lh 



(52) 



where a\h — nih/n-Q = 0.277 determined experimentally 
In the limit of infinitely high mobility of holes in the 
doped region, M\h ^> 1, Eq. |52|) goes over to Eq. (16), 



which may be written in the dimensionless form as C — 
aihd/(5 + 1). For the experimental values of the hole 
mobility, the difference between Eqs. (16) and (52) is 
about 2% . Boundary conditions in the undoped region 
are, by definition, Eq = 1, «o- = ao+ = 1, cto/i = 0. 



C. Asymptotic behavior in the doped region 

We next consider the asymptotic behavior in the spe- 
cific zones of the front: in the doped region, the undoped 
region and the transition point between doped and un- 
doped matter. We start with the doped region. Within 
the limit of high hole mobility, /j,hi 3> M± i- e - M^i ^> 1, 
we have E\ = 0, = etih, «i- = 1 + aih$/{$ + 1): 
«i + = 1— a\h/{8+l) in the doped region at £ = — oo. We 
are interested in the asymptotic approach to the doped 
state, and we investigate small deviations from the lim- 
iting values Ei = e, ah = olxh + ah, a- = ai- + cL, 
a + = ai + + a + . We also take into account the Nernst 
term (ip = 1) in the doped region, which leads to 



7e// 



-M h 



at, 



1 t 
fa 



Xh 



ah 



(53) 
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in the definition Eq. (50 1. We immediately see that, 
in the case of / = 1, Jeff diverges in the doped region 
as £ — > — oo. Thus, we have to consider two separate 
cases of / > 1 and / = 1 yielding an exponential and a 
power-law approach to the limiting values, respectively. 

We begin with the first case, / > 1, for which j e ff is a 
constant coefficient in Eqs. ( 45 1-( 48 ) , j e ff = —Mh/(f — 
1). Linearizing Eqs. ( 45 1-( 48 ) witn respect to small de- 



viations in the doped region, we obtain 



1 da h 



*lh 5 - 



da- 



da+ 



f-l d£ 8 + 1 



aihS 



aih.5 . 



5 + 1 J 8 + 1 



= — d+8 



1 - 



ot\hO _ 



8 + 1 J 5+1 



d+ — d- + (3/j. = 0. 



(54) 



(55) 



(56) 



(57) 



The system of Eqs. ( 54 )-( 57 ) has an exponential solution 
in the form e cx dh oc a_ oc <S + cx exp(x£) with positive 
eigenvalues x corresponding to the exponent decay at 
£ — y — oo. In general, the eigenvalue may be calculated 
numerically. Still, we can find an analytical solution to 
the system in the most important limit of / — 1 -C 1, 
which describes the most interesting features of the front 
in the doped region. Substituting the deviations of elec- 
tric field e from Eq. (54) into Eqs. (55) and (56), we 
obtain 



da- 



1 



1 



day. 



uih.8 

8 + 1 J (f-l)a lh d£ ' 



1 dctj 



8 d£ 



= — a+ 



1 - 



ot\h 



day. 



8+1 J (f-l)a lh d£ 



(58) 



(59) 



Thus, in the case of / — 1 <C 1, taking the difference of 
Eqs. (55) and (56) and accounting for Eq. (57) we find 
that 



ah 



5+1 + a lh {8 - 1) da h 
8(f -l)a lh dli 



with the eigenvalue 



X 



5aih(f-l) 
8 + l + a lh (S -1)' 



(60) 



(61) 



For the experimental data a\h = 0.277 and 8 = 2.2 we 
find that relaxation of the parameters to the saturation 
values in the p-doped region occurs on a length scale 
w 7.1L p /(f - 1). Therefore, in the limit of / - 1 < 1, 
the relaxation happens on length scales much greater 
than the characteristic length L p , Eq. (33), related to 



ion diffusion. We stress also that without employing the 



injection barrier in the doped region [i.e. taking -0 = in 
Eq. (40)], we do not find any deviation mode vanishing 
asymptotically to £ — > — oo in agreement with the pre- 
vious numerical results^. The physical meaning of this 
effect was explained above in Sec. IV A. 

The relaxation process goes even slower in the specific 
case of / = 1, when Eq. (54) becomes intrinsically non- 
linear 



1 da h _ ai h 8 _ 
&h d£ 8 + 1 



Then, instead of Eq. ( 60 ) , we obtain 



dan 
~dj 



5+l + a lh {8- 
with the asymptotic solution 



a h 



8 + l + a lh {8-l) 

Is 



(62) 



(63) 



(64) 



for £ — > — oo. According to Eq. (64), in the case of / = 1 



the relaxation of the hole concentration goes inversely 
proportional to the distance from the doping front. Such 
behavior is much slower than the exponential law pre- 
dicted by Eq. ([60} for / > 1. 

Thus, we obtain a smooth relaxation of the p-front 
parameters to the final values in the doped region on the 
length scales much larger than L p , Eq. (33), determined 
by ion diffusion. 



D. Asymptotic behavior in the undoped region 

In this subsection we study the asymptotic front struc- 
ture in the undoped region (£ — > oo) with eq = 1, ao_ = 
a a+ = 1, a h — 0, ip — in the limit j e ff = M m < 1, 
the latter due to the strongly nonlinear dependence of 
hole mobility on concentration, Eq. (43). Using again 
a tilde to denote the deviation variables, Eqs. (45)-(48) 



may be linearized with respect to small deviations to give 



hQ ' 



da h 



-an, 



da.- 
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a lh S 
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aihd _ 



<5 + l 



aihS 
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a- +a.h = 0. 



(65) 



(66) 



(67) 



(68) 



The system (65)- (68 1 has two independent modes in the 



form of e oc dh oc d- oc <S + oc exp(x£) decaying exponen- 
tially at £ — > oo with x < 0. One mode is related to the 
perturbations of the holes, dh ^ 0, with the eigenvalue 



X 



-1/Mj 



hOt 



(69) 



which is obtained in a straightforward way from Eq. ( 65 1 . 
In the limit of ultimately low hole mobility in the un- 
doped region, M^q -C 1, this mode is characterized by 
extremely sharp gradients, thus leading to a steep head 
of the front with an associated length scale sa MhoL p , 
which is much smaller than the length scale L p due to 
ion diffusion. This mode is also expected to lead to a 
sharp peak in the electric field, which is proportional to 
the gradient of hole concentration according to 



dan 



(70) 



The second mode in the undoped region happens with 
zero deviations of the hole concentration, &h = 0, and 
equal non-zero deviations for the ion concentrations, 

d_ 



&+ ^ 0. Then Eqs. <|66j), ([67] yield 
da± 



1 



-a± 



with the eigenvalue 



X = 



1 + 6 
2 



(71) 



(72) 



The second mode is controlled by ion diffusion with the 
typical length scale comparable to L p . Because of the 
second mode one should expect a non-monotonic behav- 
ior of the ion concentrations in the doped region. After 
sharp changes related to the first mode with complete 
vanishing of the hole concentration, a relatively slow re- 
laxation of ions to ao- = «o+ = 1 is thus expected. 



E. Behavior in the transition point 



The solution to the system (45 )-< 48 ) demonstrates that 
the doping front possesses one more specific point, which 
indicates transition from the doped to undoped zones. 
This point gives the answer to the question where the in- 
jection barrier (the Nernst term) should be switched off 
via the phenomenological parameter ip going over from 1 
to 0. A priori, it is natural to expect the critical point to 
be in the region where the hole mobility becomes com- 
parable to the ion mobility. However, it turns out that 
the system specifies an exact position of the critical tran- 
sition point (which we label by c) characterized by zero 
derivatives of Eqs. (|45l-(|48|), i.e. 
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FIG. 3. Concentrations of holes and ions, and the electric 
field inside the p-doping front; / = 1.1. 



Othc 



0. 



(76) 



Due to the essentially nonlinear dependence of hole mo- 
bility on concentration, see Eq. (|43|), the set of equa- 



tions (73l-(76) is also strongly nonlinear and it can only 
be solved numerically. Thus, solving Eqs. ( 73 1-( 76 ) 
for the experimentally obtained parameters a\h — 0.277 
and 8 — 2.2, we find the critical point at cthc ~ 0.092, 
Mh(ahc) = [i>h/fi>- ~ 8.5 with the scaled electric field at 
that point calculated as e 



(731 



0.65 from Eq. 

The critical point has an interesting physical mean- 
ing. Going back to the dimensional equations for the 
front structure, Eqs. (30)- (32 1, and setting all derivatives 



equal to zero, we find the critical point corresponding to 

- U p + n h E = 0. (77) 

Thus, at the critical point we obtain holes moving locally 
with the same velocity fJ-hE as the p-doping front, which 
may be also interpreted as a "resonance" between the 
light charges and the front. To the left of this point 
(in the doped region) hole mobility may provide faster 
velocity of the particles in comparison with the front. To 
the right of the critical point, hole mobility is too low for 
holes to keep up with the front velocity. Therefore the 
critical point plays the natural role of a border between 
the doped and undoped regions. 



= 1 



aihS 



a r - 1 



= o — - — - — a c+ o — e ( 



(74) 



(75) 



Numerical solution for the front structure 



In this section, we solve Eqs. (45)- (50 1 numerically 
with the boundary conditions ahead of the doping front 
and behind the front as obtained in Sec. IV B. Figure 
3 shows the internal structure of a stationary p-doping 
front in terms of the normalized concentrations of holes, 
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FIG. 4. Concentration of holes and the electric field inside 
the p-doping front for different values of the /-factor. 



(a) 



positive and negative ions and the electric field as a func- 
tion of distance. The front propagates to the right con- 
verting the undoped semiconductor material with low 
conductivity to the doped one with high conductivity. In 
agreement with the analysis of Subsections IV C-E, the 
whole doping front has a complicated nonlinear struc- 
ture with several specific zones characterized by differ- 
ent length scales. The length scale variations within the 
front are related, first of all, to dramatic changes of the 
hole mobility by three orders of magnitude from the un- 
doped to doped regions. We point out a sharp head of 
the doping front in the undoped region, an extremely 
smooth relaxation tail in the doped region and an addi- 
tional plateau connecting these two regions. The length 
scale variations inside the doping front resemble a sim- 
ilar effect encountered in laser deflagration, where the 
length scale may also change by several orders of mag- 
nitude within the deflagration front due to electron heat 
conduction increasing strongly with temperature^. 

As explained in Sec. IV C, the smooth relaxation tail 
in the doped region is due to high mobility of holes in 
that region and the Nernst term modeling the injection 
barrier. The numerical solution shown in Fig. 3 uses the 
parameter value / = 1.1, which provides the characteris- 
tic length scale 7lL p of the relaxation. Indeed, in Fig. 
3 we see that the relaxation length scale in the doped 
region exceeds the length scale of ion diffusion L p by ap- 
proximately two orders of magnitude. Figure 4 shows 
modifications of the front structure caused by changing 
of the parameter /. In agreement with the theoretical 
predictions, relaxation to the final doped state becomes 
smoother as / approaches unity. At the same time, the 
parameter / does not influence the head of the front. 

The details of the p-front structure in the undoped re- 
gion (head of the front) are shown in Fig. 5. Figure 
5(a) presents the concentration of holes and the electric 
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FIG. 5. Structure of the p-front in the undoped region, (a) 
hole concentration and the electric field, (b) ion concentra- 
tions. 



field. In agreement with the asymptotic theory of Sec. 
IV D, we can see sharp gradients of the hole concentra- 
tion with characteristic length scales noticeably smaller 
than the length scale L p of ion diffusion. The theory 
predicts that the electric field in the leading part of the 
front is proportional to the spatial derivative of the hole 



concentration, see Eq. ( 70 ) , which leads to a sharp peak 
in the electric field clearly visible in Fig. 5(a). Figure 
5(b) shows the variations of the ion concentration at the 
head of the doping front. Unlike the hole concentration, 
the concentrations of ions demonstrate a more compli- 
cated behavior. First, the two ion concentrations ap- 
proach each other on a short length scale related to the 
mode with the eigenvalue given by Eq. (69 1. However, 
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FIG. 6. Hole concentration and mobility inside the p-doping 
front obtained by integration from the doped region (left) with 
the Nernst term and from the undoped region (right) without 
the Nernst term. 



instead of a monotonic relaxation to the limiting values, 
the concentration n + overshoots no, so that n_ and n + 
meet at some value exceeding no (in agreement with the 
constraint of quasi- neutrality) . After that, both the ion 
concentrations n_ and n + approach the limiting value no 
together from above. Such a specific way of relaxation 
for the ion concentrations is in agreement with the sec- 
ond mode, as given by Eq. (72 1, predicted analytically 
in Sec. IV D. The ion relaxation is of course determined 



by ion diffusion and occurs on length scales comparable 
to L p . 

The specific behavior of the doping front parameters 
close to the critical point is another interesting feature 
of the system. The critical point was not discussed in 
Refill, since in that work a matching of the concentration 
gradients was done by a linear extrapolation from the 
doped and undoped regions. However, we stress that the 
critical point is not a mathematical artifact of the phe- 
nomenological transition function ip. As demonstrated in 
Sec. IV E, the local velocity of holes produced by their 
mobility is in resonance with the p-front velocity in the 
critical point, and this resonance explains the physical 
origin of the critical point and the zone around it. In 
order to understand the effect of this critical zone better, 
we may integrate Eqs. (45)- (50 1 in two opposite direc- 
tions: from the doped region with ip = 1 (from left to 
right in Fig. 3) and from the undoped region with ip = 
(from right to left in Fig. 3). Figure 6 presents the con- 
centration of holes and the hole mobility obtained from 
such an integration procedure. As illustrated in Fig. 6, 
integrating from the doped region we could not reach the 
undoped one, and vice versa. Instead, using both ways 
of integration we reach asymptotically the critical point 
as a saturation limit: from above for -0 = 1 and from 



below for ip = 0. From the mathematical point of view, 
passing the critical resonance point means changing sign 
of the right-hand side of Eq. (EBJ, (CM h e-l). The 
continuous transition from the doped (left) to undoped 
(right) regions requires a non-positive derivative of the 
hole concentration, dah/dt;, and, therefore, the effective 
hole diffusion 7 e // has to change sign in the critical point. 
We remind that the effective diffusion is a combination 
of the real diffusion and the Nernst contribution describ- 
ing the injection barrier. The role of the Nernst term 
is controlled by the phenomcnological function ip, which 
therefore has to change from 1 (the Nernst term is then 
switched on) to (the Nernst term is switched off) in the 
vicinity of the critical point. For example, matching the 
concentration profiles obtained by integration from the 
right and from the left at the critical point we find the 
front structure for ip(rih) in the form of a step- function. 
Still, a smooth transition function is required for treating 
the doping fronts numerically within the evolution prob- 
lems. Due to the lack of a good thermodynamic model 
for the transition from the undoped to doped state, we 
use the following phenomenological form of the transition 
function: 



ip = 0.5 + 0.5tanh [A x (CM h e - 1) \CM h e - 1| 
+0.51n(a/, oo /a/ i - 1)]. 
The hyperbolic tangent of Eq. (78) provides a smooth 



(78) 



transition from 1 to for ip as we go from the doped to 
undoped region. The function Eq. ( |78[ ) depends on the 
combination (C^e — 1), which is the right side of Eq. 
( |45| . When this combination is zero, the residual term 
makes 7e// = and allows for a smooth change in all 
the concentrations. The other parameters of the transi- 
tion function are chosen to reduce the plateau near the 
critical point. In our calculations presented in Fig. 3 
we used A\ = 8 and n = 0.1. Hypothetically, a transi- 
tion function ip may exist, which eliminates the plateau 
completely. However, we believe that it is important 
to obtain the transition function from first principles of 
thermodynamics and quantum mechanics, rather than to 
make a more elaborate phenomcnological construction. 
This is indeed an important and difficult problem left 
for the future, while at present we simply use a phe- 
nomenological function for ip to obtain a smooth transi- 
tion through the critical point from the undoped to doped 
region. 

All the main characteristic features of the doping front 
for holes are relevant for electron doping as well. In Fig. 
7 we depict the stationary doping front for electrons (the 
n- front). The characteristic length scale related to the n- 
doping front is defined in the same way as for the p-front 
according to 

L n = — = D ~ Ule (79) 
U n0 (p+ + fi-)n E Q ' 

In order to obtain a numerical result for the front struc- 
ture we used the electron mobility function 

M h = 145 [1 + tanh (21.6a/, - 4.3)] (80) 
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FIG. 7. Concentrations of electrons and ions, and the electric 
field inside the n-doping front. 



FIG. 8. Concentrations of holes and ions, and the electric 
field inside the p-front of weak doping for n lh /no = 0.08. 



constructed from the experimental dat a 113 * 30 -*, and the 
phenomenological transition function similar to the one 
given in Eq. 



(78) 



ij; e = 0.5 + 0.5tanh 



14.4 (CM e e- 1)| CM e e- l| c 



+0.5 In (a eoo /a e - 1)]. 



(81) 



The electron mobility in OSCs is somewhat lower than 
the hole mobility. Similar to the p-front shown in Fig. 3, 
the n-doping fronts have also very elongated tails in the 
doped region, a sharp head with a peak of the electric 
field in the undoped region and a critical point of tran- 
sition from the undoped to doped state at /i e //i_ w 3.5. 
We point out that the critical point for electrons corre- 
sponds to considerably lower mobility (the critical point 
for holes is achieved at (J,h/v>- ~ 8.5). Because of the 
lower electron mobility, the critical point demonstrates 
a more complicated structure for the n-front in compar- 
ison with the p-front. In the case of the p-front, satu- 
ration of all the concentrations and the electric field to 
the plateau of the critical point occurs monotonically, 
see Fig. 6. On the contrary, in the case of the n-doping 
front presented in Fig. 7, the concentration of positive 
ions and the electric field exhibit a non-monotonic behav- 
ior when approaching the plateau. This behavior shows 
clearly existence of two perturbation modes in the doped 
zone close to the critical point, which resemble qualita- 
tively the modes obtained in Sec. IV D. We have found a 
qualitatively similar structure of the transition zone for 
several types of function if). 

An interesting consequence of the above discussion of 
the critical point is the possibility of a weak doping pro- 
cess, when the final concentration of the light charge car- 
rier is still smaller than the corresponding critical value. 
Such a doping front can be described without the Nernst 
potential and the transitional function. At the same 



time, the front structure in the undoped region remains 
the same as described in Sec. IV D. The p-front of such a 
weak doping transformation is shown in Fig. 8; it is sim- 
ilar to the right part (head) of the front in Figs. 3 and 6. 
The final hole concentration in Fig. 8 is smaller than the 
critical one, namely, a\h — 0.08, and the corresponding 
final mobility of the holes is fiih/f^o = 4.5. As the hole 
mobility is about 80 times smaller than in fully doped 
case, then, according Eq. ( 77 1 , the final electric field re- 
mains relatively large behind the front, e\ ~ 0.8, though 
smaller than the initial one, eo = 1. The total length of 
the front is much shorter than the front width in Fig. 3, 
as the long relaxation tail in the strongly doped region is 
missing here. In general, the front of weak doping trans- 
formation may be interpreted as a part of the complete 
doping front from the head to the critical point. 

Finally, we present the structure of the p- and n- 
doping fronts as they accelerate towards each other in 
Fig. 9. Since the characteristic time scales related to 
the fronts (that is r p cx L p /U p = D_/U p and t„ cx 
L n /U n = D-/Un) are much smaller than the time of the 
front acceleration, then the structure may be obtained 
within the quasi-classical Wentzel-Kramers-Brillouin ap- 
proximation. Within this approximation, structure of the 
fronts remains self-similar, but the length scales L p and 
L n decrease due to the increase of the electric field and 
the front velocities U p and U n as discussed in Sec. III. 



V. CONCLUSIONS 

In this paper we have investigated the dynamics and 
internal structure of the planar p- and n-type doping 
fronts in organic semiconductors with applications to 
LECs. The study is based on the drift-diffusion model 
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FIG. 9. The p- and n- fronts accelerating towards each other, 
(a) hole and electron concentrations are depicted by magenta 
and green, (b) negative and positive ions are shown by blue 
and red. 



of processes/devices. In the doping front (studied in the 
present paper) holes and electrons populate the active 
material while ions give way to the light charges. As a 
result, conductivity of the material increases drastically, 
by 2-3 orders of magnitude. On the contrary, the trans- 
formation fronts in ambipolar devices imply that cathions 
replace holes with strong decrease of conductivity. The 
different characters of the processes lead naturally to dif- 
ferent properties of the transformation fronts. 

Here we have studied parameters of the doping fronts 
on different scales, both within the discontinuous front 
approach and also by taking into account the internal 
front structure. Within the limit of a discontinuous front 
we have derived the analytical formulas for the front ve- 
locities, the ion concentrations in the doped region, and 
described dynamics of the p- and n-doping fronts acceler- 
ating towards each other in LECs. The analytical results 
for the planar front dynamics are in a good quantita- 
tive agreement with the experimental data for the slow- 
est part of the experimentally observed front brush. We 
remind that experiments demonstrate also a complicated 
multidimensional front dynamics related to instabilities. 
Theoretical investigation of the front instabilities is be- 
yond the scope of this paper and is presented elsewhere^. 

One of the main purposes of the present paper was 
to investigate the internal structure of the doping fronts. 
In agreement with the previous ideas^, we show that 
continuous transition from the doped to undoped state 
in the form of a moving front requires a thermody- 
namic injection barrier and a nonlinear dependence of the 
hole/electron mobility on concentration, which is quite in 
line with the basic principles of the doping process. We 
have studied the asymptotic behavior of the front param- 
eters: 1) relaxation to the doped state at the back of the 
front, 2) deviation from the undoped state at the head of 
the front and 3) the critical point of the transition from 
the doped to undoped parts of the front, where the ve- 
locity of the light charges are in resonance with the front 
speed. We have also obtained a numerical solution for 
the front structure. Both the analytical theory and the 
numerical solution demonstrated the multi-scale features 
of the doping fronts, which include an extremely smooth 
relaxation tail in the doped region, a sharp head of the 
front with large gradients in the undoped zone, and a 
plateau at the critical point. The described front struc- 
ture agrees qualitatively with observations of the previ- 
ous experiments^. 



taking into account the injection barrier and a strongly 
nonlinear dependence of the hole and electron mobil- 
ity/diffusion on concentration, following the work Ref. 
l28l A similar model has been employed before to de- 
scribe transformation fronts in one-electrode devices like 
electrochemical sensors and actuators in Ref. [37J Still, 
there is an important difference between these two types 
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